The Heating of Thermal Electrons in Fast Collisionless Shocks: 
The Integral Role of Cosmic Rays 

Cara E. Rakowski^ 

Space Science Division, Naval Research Laboratory Code 7674R, Washington, D.C. 20375 

O . 

O ■ J- Martin Laming 

(N 

Space Science Division, Naval Research Laboratory Code 767 J^L, Washington, D.C. 20375 

a _ 

^ ^ and 

^ ■ Parviz Ghavamian 

' (-('' Department of Physics and Astronomy, Johns Hopkins University, 34 00 N. Charles Street, 

P^- Baltimore, MD, 21218 

O 

parviz@pha.jhu.edu 

C/3 ' 



> 

00 



ABSTRACT 



Understanding the heating of electrons to quasi-thermal energies at collision- 



O ' less shocks has broad implications for plasma astrophysics. It directly impacts 
CO ' 

^ ■ the interpretation of X-ray spectra from shocks, is important for understanding 

Q ■ how energy is partitioned between the thermal and cos mic ray populations, and 

oo ■ 
O 



provides insight into the structure of the shock itself. In lGhavamian et al.l (120071 ) 
we presented observational evidence for an inverse square relation between the 
■ electron-to-proton temperature ratio and the shock speed at the outer blast waves 

;_i ' of supernova remnants in partially neutral interstellar gas. There we outlined how 

^ lower hybrid waves generated in the cosmic ray precursor could produce such a 

relationship by heating the electrons to a common temperature independent of 
both shock speed and the strength of the ambient magnetic field. Here we ex- 
plore the mechanism of lower hybrid wave heating of electrons in more detail. 
Specifically we examine the growth rate of the lower hybrid waves for both the 
kinetic (resonant) and reactive cases. We find that only the kinetic case exhibits 
a growing mode. At low Alfven Mach numbers (~15) the growth of lower hybrid 
waves can be faster than the magnetic field amplification by modified Alfven 
waves. 



^National Research Council Fellow 
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Subject headings: ISM: Cosmic Rays, ISM: Shock Waves, ISM: Supernova Rem- 
nants 



Introduction 



The main accelerators of cosmic rays (CRs) are widely believed to be high Mach number 
shocks in coUisionless plasma, here loosely defined as plasma where charged particles interact 
predominantly through plasma wav es rather than by Coulomb collisions (for a thorough 
conceptual and historical review see iMalkov fc O'C Drurvl[2nnil ). However, a concensus is 



emerging that CRs are not simply a byproduct of coUisionless shocks, but in fact play an 
integral role in the shock structure, dynamics, and energetics. For example sound waves in a 
cosmic ray pressure g radient can smooth out the shock jump in cosmic ray modified shocks 



(jPrury fc Falldll986l ). More recent analytic and numerical work has shown that modified 



Alfven waves in the CR precursor may amplify the magnetic field to many times its ambient 



Lucek & Bell 


( 


2000 


); 


Bell & Lucek 


( 


2001 


); 


Bell 


( 


2004. 


2005) 



Observational support for 
dramatic magnetic field amplification ahead of shocks exists in the form o f extremely thin 



X-ray synchrotro n rims of supernova remnants (SNRs) su ch as Cassiopeia A (jVink fc Laming 



2003 


), SN 1006 ( 


Lone; et al. 


2003: 



20041. and Tvcho's SNR (Warren et al 



2005; Cassam-Chenai et al. 2007). As noted by ICassam-Chenai' et al.l (120071 ) there are two 
possible interpretations for the narrow width, however both require a dramatically amplified 
magnetic field ahead of the shock. Either the rims are thin because the high magnetic field 
causes rapid synchrotron cooling of the X-ray emitting electrons or the scale of the rims 
represents the scale of magnetic field de-amplification behind the shock. 

In this paper we explore another area where CRs may infiuence the properties of the 
shock, namely through the heating of quasi-thermal electrons. For shocks in coUisionless 
plasma the heating of electrons must occur through the damping of waves generated by the 
other more massive charged particles that dominate the energetics. Given the wide array of 
possible plasma instabilities at coUisionless shocks an observational relationship for electron 
temperature, Tg at the shock front was required to limit theoretical discussions. An inverse 
relationship between the initial ratio of electron to proton temperatures immediately behind 
the shock, (T^/Tpjo , and the shock velo city, v,^ , has been reported in a series of observational 



pape r s on supernova remnant shocks (IGhavamian et al.l 1200 ll . 120021 . l2003l : iRakowski et al. 



20031 : iRakowskil l2005l : IGhavamian et al.l |2007| ) and has also beei i noted among the h igher 
Alfven Mach number e vents in a sample of solar wind shocks (ISchwartz et al.lll988l ). In 



Ghavamian et al.l (120071 ). we focused on shocks propagating into partially neutral gas. Here 
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the collisional excitation of broad and narrow Balme r line e mission at the shock fron t can 
be used to diagn o se (T p/Tp)n- In iGhavamian et al.l (I2OOII ). iGhavamian et ahl (120021 ) and 
Ghavamian et al.l (120071 ) we described the method of simultaneously constraining Vg and 
(Te/Tp)o via measurement of the width of the broad Balmer line and the ratio of broad 



to narrow Balmer line flux (see also iHeng fc McCrayl 120071 : iHeng et al.l 120071 ) . Our results 
are consistent with an inverse square relationship, (Te/Tp)ooc l/v^ for shock speeds above 
~400 km s~^ (Ghavamian, Laming & Rakowski 2007). Given that Tp oc at the shock front 
by the Rankine-Hugoniot jump conditions, the inverse relationship between equilibration 
and shock speed implies that the electron temperature itself is nearly constant ~0.3 keV, 
independent of shock speed. 

The insensitivity of electron temperature to shock velocity suggests a heating mechanism 
within the extended diffusive cosmic ray precursor ahead of the shock. In this case the 
electron heating would be more reflective of the generic properties of cosmic ray acceleration 
and diffusion than tied to the specific attributes of the shock. In contrast, prior work 
on heating by shock- reflected ions that are confine d to within a gyroradius of the (quasi- 
perpendicular) shock (ICargill fc PapadopouloslllQSSl ) sugge sted that (Te/Tp)o would remain 
constant with shock velocity. In IGhavamian et al.l (120071 ) we suggested that lower hybrid 
waves in the cosmic ray precursor of a perpendicular shock might be a plausible electron 
heating mechanism. 

Lower hybrid waves are electrostatic ion waves directed nearly perpendicular to the 
magnetic field with a frequency equal to the geometric mean of the electron and ion gyrofre- 
quencies. Electrons that would otherwise screen the ion oscillation are pinned to magnetic 
field lines. In addition, the group velocity parallel to the magnetic field greatly exceeds the 
group velocity perpendicular to the field (cu/A^ji >> uj/k±). Therefore, the wave can simul- 
taneously resonate with ions moving across the field lines and electrons moving along the 
field lines, facilitating coUisionless energy exchange between them. Based on simple argu- 
ments about the width of the cosmic ray precursor and the electron diffusion along the field 
lines, we showed that electron heating from lower hybrid waves in the cosmic ray precursor 
would be independent of both the shock speed and the magnetic field. Here we explore this 
mechanism in more physical detail. 

In Section 2 we calculate the growth rate of lower hybrid waves, first examining the 
kinetic (resonant) then the reactive (non-resonant) case. The treatment h ere is mathemat- 
ically similar to the work on modified Alfven waves by lAchterberg] (119831 ) and iBelll (12004 . 
2OO5I ) involving the cosmi c ray contribut ion t o the plasma die lectric tensor. The analysis 
also draws on the work of iLamingl (l2001al ) and iLaming (l2001bl ) on lower hybrid waves from 
shock-refiected ions. We compare these growth rates with those for magnetic field amplifi- 
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cation, to assess the conditions under which electron heating might occur. In section 3 we 
discuss the structure of the cosmic ray shock precursor in more detail. We pay particular 
attention to the magnetic field geometry, since the excitation of lower hybrid waves requires 
a quasi-perpendicular shock. We show schematically how magnetic field amplification and 
lower-hybrid wave heating might co-exist in the shock precursor for either parallel or per- 
pendicular initial geometries. We also review some other ideas for electron heating, and 
make some quantitative predictions from our model for various shock parameters. Included 
in the appendices are a discussion of cosmic ray diffusion coefficients and a derivation of the 
resonant growth rate for electromagnetic waves. 



2. Cosmic Ray Growth Rate of Lower Hybrid Waves 



Lower hybrid waves ahead of collisionless shocks have particularly interesting properties. 
They can have a group velocity away from the shock equal to the shock velocity itself 
(IMcClements et al.lll997l ). This can in principle allow the waves to grow to large amplitudes, 
even if their intrinsic growth rate is small. To determine if lower hybrid waves can heat the 
electrons to the ~0.3 keV temperature observed, we must first calculate the growth rate of 
this instability to see if it will have sufficient power to overcome the damping effect of the 
electrons as well as to compete with other instabilities in the precursor. We calculate this 
growth rate in both kinetic and reactive limits, i.e. either considering the cosmic rays with 
energies in resonance with the lower hybrid w ave frequency or the integrated contribution 
of the entire distributio n, respectively ( see e . g iMelrosd Il986l) . Related kinetic and reactive 



calculated in iLaming J2001ah and lLaminel J2001bh . but only for the case of shock 



cases were 

reflected, non-relativistic ions gyrating around the magnetic fleld, represented as a particle 
beam. Here we begin the discussion with the resonant case. 



2.1. Kinetic Growth Rate 



We model the normalized cosmic ray distribution function diffusing upstream as 



fip) 



n 



CR 



(2K-3)r(K) 



4v^(7r/€)3/2p3 r(/€-l/2) 



1 + 



{p^ - mvs)^ + pI + pI 



2kp- 



where the coordinate system is aligned so that the shock speed Vg lies in the x-direction. 
Here pt is deflned as the "thermal" momentum, which we take to be (3/4) mf^, and n[ 



CR 



denotes the density of suprathermal particles with distribution function / (p). The functional 
form above, known as a "kappa" distribution, is often seen for particle distribution functions 
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associated with shocks in the solar wind, and may be derived as e quihbrium distributions for 



a system of particles and waves under certain conditions (see e.g. iLaming fc Lepril 120071 ). in 
contrast to a system of particles only which gives a Maxwellian distribution of width p^. The 
kappa distribution resembles a Maxwellian for p < pt, and in fact for k ^ oo is exactly a 
Maxwellian. At higher particle momenta it tends smoothly to a distribution / oc p"^'^. Below 
we shall take k = 2 to model the well known f{p) oc cosmic ray distribution predicted 
by diffusive shock acceleration in shocks with a compression ratio of 4. In connecting the 
cosmic rays to the lower energy particles in this way, we are somewhat blurring the distinction 
between "cosmic rays" and other suprathermal particles reflected from the shock. Hence we 
denote the combined density of these particles as nj^^ to distinguish it from density of true 
cosmic rays, ucr, which will appear in expressions derived by other authors. Note that 
all these particles are distinct from the ambient thermal plasma upstream of the shock, 
which here is considered to be a Maxwellian with much lower temperature than pt used in 
Equation (1) for the upstream suprathermals. When discussing the kinetic instability we 
will be focusing on particles that obey a diffusion equation ahead of the shock, due to their 
interaction with turbulence, rather than gyrating around field lines. We will qualify this 
distinction more carefully below in our discussion of the reactive instability. 

The appropriate dispersion relation for Equation (1) can be found from the cold plasma 
dieletric tensor. For electrostat ic waves at frequencies close to the lower hybrid wave fre- 
quency we have (Laming 2001al jbl) 



where 6 denotes the angle between the wavenumber of the perturbation and the preshock 
magnetic field; = {Anq^nx / rrix)'^^/'^^ is the plasma frequency of a given species x (electrons, 
ions, cosmic rays...) with charge g, density n^;, and mass m^; VL^ = qB / (m'jLc) is the 
cyclotron frequency of species x (with being the Lorentz factor); and unadorned u being 
the lower hybrid wave frequency which is the geometric mean of the electron and proton 
cyclotron frequencies. Using the Landau prescription for evaluating the integral at the 
resonant pole and taking only the imaginary parts of the dielectric tensor equation, we find 
the growth rate for the lower hybrid waves 



To compute 7, we take n'(j^ oc exp —xvg/ D where D is the cosmic ray diffusion coef- 
ficient and is in principle dependent on the cosmic ray momentum (making / = D/vg the 
characteristic diffusive lengthscale) . With this substitution we begin the evaluation of the 
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last integral in Equation (3) 



{2k-?,)T{k) 2 



TXK 



1 + 



{p^ - mvsf 



(4) 



4v^(7r/c)3/2p3 V{k-1/2) K-l 

where we have separated out the components of p perpendicular to the shock (and k) from 
Px- Substituting back into Equation (3) yields a growth rate 



7 



7r\3/2 q 



{2k - 3)T{k) 



X < - 



1 + 



'2 +Cj2^COs2e V2r(/t- 1/2) 

{px - mvs)^] '"^ {pa: - mf. 



S{uj — k • v) 



2 up 



Pt 



+ 



K — 1 D"^ dpx 



2kp 



1-k" 



xVs/D 



For waves to stay in contact with the shock w/k c::^ —2vs f lLamingll2001al ) in the cold plasma 
electrostatic limit, i.e. reflected ions returning to the sh ock excite the w aves. This remains 
generally true when these approximations are relaxed (ILamingl l2001bl ). so the 6 function 
picks out Px = —2mVs (cosmic rays returning to the shock). Hence 
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7r 
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2kp- 
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^/2V{k - 1/2) 
KxVs dD 

{K-l)D^d^x 



~2^ 
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(6) 



Assuming dD/dpx = (see Appendix A on cosmic ray diffusion coefficients) and k = 2, only 
the first term within curly brackets remains. Evaluating the F functions and substituting 
Pt = (3/4)mfs we arrive at the following expression for the kinetic growth rate of lower 
hybrid waves 

8 

7 



225 iul + col cos2 e 



-xvs/D 



(7) 



pi 



pe 



where 0;^ denotes the plasma frequency for nj^^. Substituting in the frequency definitions 
we note that approximately 7 oc (n^^/ni)f2i. Before proceeding, we pause to compare this 
growth rate with those for magnetic fi eld amplification. In the case where Alfven waves are 
resonantly excited, the growth rate is (jMelroselll986l : iPelletier et al.ll2006l ) 



7s,i 



3n fli ricR f cos 9 
16 va ni \ \ cos 6 



■Vs COS ( 



^va - '^Vs sin0 ]k\\rg 



(8) 



where va is the Alfven speed, rii is the density of ions in the background plasma and is 
the gyroradius of cosmic rays. This expression differs from that in the cited references in the 
factor cos (f) and the term in sin 0, where is the angle between the shock velocity and the 



- 7- 



magnetic field. At perpendicular shocks, the growth rate of resonant Alfven waves can be 
neglected, but at parallel shocks may be larger than that for lower hybrid waves, depending 
on the ratio n'(j^/ncR. However as we shall argue later, all shocks subject to magnetic field 
amplification become perpendicular, and this is the geometry whe re lower hyb rid waves are 
most effectively excited, so we neglect '~^B,res from here onward. iBelll (120041 ) discovered a 
nonresonant growth rate for Alfven waves, with approximate growth rate 



IB; 



ncR 



Hi 



(9) 



which has a maximum value of MAQincR/2ni. According to iBelll (120051 ). this instability 
operates for arbitrary orientations of B, Vg and k, indicating that it will also amplify magnetic 
field at perpendicular shocks. Its growth rate is strongest for k||B and zero for k ± B. 
Equating the maximum value of 'jB,nonres with the lower hybrid wave growth rate calculated 
above, we find the critical Alfvenic Mach number ~ Sn'^jj^/ricR, such that for higher 
Myi, cosmic rays preferentially amplify magnetic field, and for lower they generate lower 
hybrid waves. The numerical value depends on the ratio n'(j^/ncR- In the next section we 
will argue that these two densities should not be the same, and that n'(jp^ > hcr, following 
from a consideration of the reactive growth rate for lower hybrid waves. 



2.2. Reactive Growth Rate 



The reactive case involves the integrated contribution to the growth rate from the entire 
cosmic ray distribution. Thus we will examine successive orders in an expansion of f{p) to 
see if they produce any growing modes. Although no instability is found in this process, we 
do uncover potentially interesting constraints on the properties of f{p). 

We consider agai n the last ter m in Equation (2), the cosmic ray contribution to the 
dielectric tensor (e.g. 
resonant case; 



Melrosd Il986l ) which includes a factor that reduced to unity for the 



K 



CR 



A-nq^ f k - df /dpk-v 



(10) 



k"^ J LU — k- V oj 

For the case of a beam of cosmic ray part icles localized around Vg, one recovers the usual 
beam reactive instability ( jLamingl l2001bl . equation (A6)). However, as we demonstrate 
below, for a more physical quasi-isotropic cosmic ray distribution no instability is recovered. 
We consider the case where cosmic rays drift with velocity Vg, and waves are generated with 
k||vs. The CR distribution function from appendix A, expanded in terms of the cosine of 
the angle between k and v or k and Vg, cos a = /i, is 



/ 



/o(l + /i^ + 

V 
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such that J fjLvfdfl — foVg as before. Neglecting terms of order Vg/v"^, the CR contribution 
to the dielectric tensor becomes 



2% 



-1 ^ 



djji—p dp. 

UJ — KVH op 



(12) 



Expanding the numerator into terms divisible by {ui/kv — ix) and evaluating the integral over 
/i we obtain 
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(13) 

We then evaluate this integral in the two limiting cases away from the pole, uj kv and 
UJ <^ kv. For UJ ^ kv; 

UJ — kv 



In 



leading to 



2^2 /-oo 



UJ 



~2vs 
kv 
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All terms in brackets cancel save for one, giving 

g^2^2 /-ooo^, 1fi^2^2 



K 
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(15) 



(16) 



Iq 3uj dp^ ^ 3uj'^ 

the first term goes to zero so long as /o(oo) faster than and the second term is 
oc ncR/i'jm) leaving 

« -^pW^'- (17) 

and 



UJ, 



11? 



uj' 1 + T^sin^ e - ut, - ujLcos' 9 - u, 



pi 



pe 



'^pCR 



0. 



(18) 



This simply modifies the l/uj"^ term in the dielectric tensor, changing the frequency of the 
solution but not creating any complex roots, hence no instability is generated. 



Likewise, in the case where uj <^ kv; 

UJ — kv 



In 



UJ + kv 



2uj 
kv 



(19) 



leading to 
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dp 
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AU but the first term are negligible in this limit, hence 

cok Jq V op 
This leads to a full dispersion relation that can be written as 

which also lacks complex roots, regardless of the actual evaluation of the integral. 

Higher order terms in the expansion of / = /o (l + {3vs/v) + fi^ {3vs/v f /2 + fi^ (3fs/f )^ /6 + ...) 
give rise to higher order terms in uu. For uu » kv, the dispersion relation equation 18 be- 
comes to next highest order 

(l + ^ sin^ 0^ - (a;^, + < cos^ 9 + u^l^n " 'i^^Vv^ J I'^p'dp) u 

+f7rVKVJfp> = 0, (23) 

which is stable since the terms oc / ^^^^p^dp are of order ~ ^pCBPlh"^ '^pCR- When 
uj «kv the dispersion relation Equation (22) takes on the next highest order terms 

- -^T- J -.^P dp (3- + - kv:) - J -,^P dp [u: + -.kv.j (24) 

to become a quartic equation. This has four real solutions so long as J ^^p'^dp > 0. In 
fact if J ■:^^p'^dp < 0, the addition of higher order terms in the expansion of the cosmic 
ray distribution function would dramatically alter the character of the solutions, a situation 
that must be considered unphysical. We require / ^^p'^dp > 0, which means that at low 
momenta, dfo/dp > 0, and the cosmic ray distribution cannot be monatonically decreasing 
from V = p = 0. Our use of the kappa distribution in the preceding section may therefore 
be questioned. However the resonance at Px = —2mVs places it well into the region of 
the distribution where the gradient is negative, and so modifications to the low momentum 
behavior would have very little effect on our result. However this observation does imply 
that the distribution of particles obeying a diffusion equation ahead of the shock is unlikely 
to extend down to zero momentum. Some natural break must exist between the quasi- 
thermal population gyrating around field lines and the cosmic rays diffusing in turbulence. 
The forgoing also neglects the cosmic ray induced current in the background plasma. The 
inclusion of such effects leads to the modification tu^j — >• 0;^^ — utpCR, and has no effect on 
reactive instabihties. 



(21) 
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Discussion 



3.1. Electron Heating or Magnetic Field Amplification? 

We have calculated the growth rate for waves that damp by heating electrons, in a cos- 
mic ra^^_^ock precursor using similar approximations and techniques to those employed by 
Belli (120041 ) ■ Both the lower hybrid wave heating of electrons and the growth of magnetic field 
through modified Alfven waves redistribute energy within the cosmic ray precursor. An im- 
portant question is which of these is more effective, i.e. which grows faster? Above we derived 
a critical Alfven Mach number. Ma — 3n'fj^/ncR, which divides the regime of magnetic field 
amplification from that of lower hybrid wave growth. Following from the treatment of the re- 
active instability above, we estimate ucr = /J^ fcR^T^p^dp ~ (S/vr) n'^j^Vs/vinj where fcR is 
given by equation 1 and Pinj is the injection momentum where particles may begin to partic- 



ipate in a diffusive shock acceleration process. The approximate result n'^R/ncR 
gives Ma — 1.5vinj/vs as the critical Alfven Mach number. 



/2vs 



The next step in determining the critical Alfven Mach number is to find an appropriate 



for the injection of seed particles into the cosmic ray acceleration process. IZank et al. 



(120061 ) argue that quasi-perpendicular shocks have similar injection requirements to quasi- 
parallel shocks, but that highly perpendicular shocks require much higher injection energies. 
In the case of nonresonant magnetic field generation, we consider the case of a highly per- 
pendicular shock since the gene rated magnetic fie ld will be perpendicular and much stronger 
than the initial magnetic field. IZank et al.l (120061 ) give the injection velocity as 

1/2 



-'in J 



3^. + ^25) 

_(r-l) _ 

where r is the shock compression ratio, and DBohm and Z)_l are the cosmic ray diffusion 



coefficients in the Bohm limit and in the perpendicular direction respectively. iReville et al. 



(120081 ) give DBohm/ D± ~ 3 for cosmic rays where kvg 



and so f j„j 



lOfo. Thus the 



Alfven Mach number at which lower hybrid wave growth takes over from magnetic field 
amplification should be about 15, unless the magnetic field saturates at a lower value (i.e. 
higher Ma) before this is reached. 

The growth of lower hybrid waves is most efficient at a quasi-perpendicular shock, 
whereas the gr owth of ma gnetic field through modified Alfven waves is strongest at a quasi- 
parallel shock (IBelll l2005l ). This a ppare n t contradiction is a ctuall y easily resolved. At an 
initially quasi-parallel shock, iBelll (l2005l ). iReville et al.l (l2008l ) and IZirakashvili et al.l (l2008l ) 
show that a highly helical magnetic field develops. The distortion of an initially parallel 
field line is shown schematically in Figure 1, showing the evolution of the shock from quasi- 
parallel to quasi-perpendicular. A similar schematic in Figure 2 shows the evolution of an 
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initially quasi-perpendicular shock, where magnetic field is amplified orthogonal to the pre- 
existing magnetic field, but where the shock remains quasi-perpendicular. In both cases a 
perpendicular field is generated, thus allowing lower hybrid wave growth and electron heating 
in a region close to the shock as indicated. 

Another poten tial proble m is the cavities seen in simulations of the growth of modified 
Alfven waves (e.g. iBelll 120051 ). The helical field from an initially quasi-parallel geometry 
naturally creates a filamentary structure, dragging the thermal plasma with it, while cosmic 
rays tend to accumulate in the low density cavities. This is problematic for our mechanism 
that requires spatial coincidence between cosmic rays, magnetic field and thermal plasma. 
A possible solution is that the growth of lower hybrid waves takes over from the growth 
of modified Alfven w aves, so th a t the cosmic r ay dr iven magnetic field nev e r rea ches its 
final saturated state. Eel (jSQol), iReville et aP J2008h and Izirakashvih et al.l koO^ derive 
a saturation magnetic field by setting 'yB,nonres = in equation (9) to give 6B ~ jrg/in 
or 6B'^/8tt ~ ncRTniVsVinj/2. This gives an Alfven Mach number at saturation of M\ ~ 
rii/lOricR, (assuming Vmj ~ lOf^), which for likely parameters rii/ricR ~ 10^ gives a value 
of Ma of similar magnitude but possibly lower than that where the electron heating is 
expected to take over. Bearing in mind that we took the strongest growth rate for magnetic 
field amplification to estimate where electron heating takes over, it is quite plausible that the 
amplified magnetic field never reaches saturation. Also, as the initial shock state becomes 
more quasi-perpendicular, the growth rate slows down, and the circularly polarized Alfven 
waves become elliptically polarized, ultimately becoming linearly polarized in the limit of a 
true perpendicular shock, eliminating the growth of such cavities. 



Pelletier et all fl2006[ ) find that the nonresonant instabihty of [Bell (l2004f ) andlBeUl (120051 1 
dominates over the more familiar resonant instability when the shock velocity Vg is greater 
than a few times ecRC where enn is the ratio of the cosmic ray energy density to the kinetic 
energy density of the shock. iNiemiec et al.l (120081 ) simulated the cosmic ray driven ampli- 
fication of magnetic field in a parallel shock using Particle-In-Cell simulations, which can 
naturally account for the backreaction of the generated magnetic field on the cosmic ray 
current. In conditions where the nonresonant mode should grow, they find magnetic field 
amplification only to 5B ~ B. The magnetic field again produces filaments, but they do not 
find cosmic ray accumulation in the filament cavities. They do not find strong growth, and 
argue that saturation occurs because the incoming fiow to the shock is decelerated by the 
cosmic rays, reducing their relative velocity and hence the cosmic ray current. 

For our electron heating model, the precise degree of magnetic field amplification is 
unimportant so long as the cosmic ray diffusion coefficient remains proportional to 1/B. It 
is only necessary that the shock be sufficiently quasi-perpendicular to allow cosmic rays to 
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generate lower hybrid waves. A reduced cosmic ray current does not necessarily produce an 
appreciable affect on the kinetic growth rate for lower hybrid waves. So long as the current 
does not vanish, the initial effect of reducing Vg in equation 6 is to bring more cosmic rays 
into resonance with the lower hybrid waves. Another estimate of the cosmic ray density 
necessary to heat electrons may come from the long wavelength limit of the magnetic field 
amplificati on, when 7^ = -^/k ■ BncRqVs/nimi, for both parallel and perpendicular cases 
( jBelll I2OO5I ). Electron heating then requires 7 ~ VtiricR/ni > ^JtlcrVs cos (p/uiVinj^^i, taking 
k = Qi/vinj (probably an overestimate), yielding ncn/ni > cos 0/10. At cos0 < ^J'mjmi^ 
the values typical for lower hybrid wave propagation, the value for ncRlrii is low enough 
(0.001 - 0.01) to make electron heating by cosmic rays plausible. 



3.2. Other Electron Heating Mechanisms 



Several other researchers ha ve considered the g e nerat ion of waves in a shock pr e cursor 
as a means of heating electrons. lOhira fc Takahara J2007h an d IShimada fc Hoshind (120001 ) 
have both considered the model of 



Cargill k. Papadopoulos Jl988h in n aore detail, using 



Parti c le- In- Cell codes ra t her t h an a hybrid approa ch. Other references (jPieckmann et al. 



2OOOI : iMcClements et al.l I2OOII : ISchmitz et al.l |2002| ) focus more on the electron injection 
problem for diffusive shock acceleration, rather than the thermal electron temperature, again 
invoking various wave modes in a reflected ion precursor. Our principle departure from these 
works has been to treat similar wave modes upstream of the shock, but excited by cosmic 
rays undergoing diffusive shock acceleration rather than by quasi-thermal ions reflected from 
the shock. This allows electron heating to occur over a much more extended upstream 
region dictated by the cosmic ray diffusion coefficient, rather than the ion gyroradius. 
In addition, expressing the thickness of this region as / ~ D/vg naturally results in elec- 
tron heating that is essentia lly independent of shock speed, as argued from observations of 
Balmer-dominated shocks in lGhavamian et al.l (120071 ). On the other hand if the thickness of 
the electron heating region is comparable to the ion gyroradius then I oc Vg and Te oc vl- 
This results in constant Te/Tj with Vg-, contrary to what is observed. 

A number of other authors have investigated the role of the cross-shock potential in 
heating the electrons. Inside the (quasi-perpendicular) shock ramp, the magnetic field may 
"overshoot" , i.e. increase to a value greatly in excess of its asymptotic downstream strength 
before decreasing again. The electric field arising from the small charge separation associated 
with this magnetic field gradient, £^ ~ 9/9a; {B^) / {Sireni), can decelerate ions and accelerate 
electrons. Such effects are know n to be important at low Mach number shocks where a 
laminar approximation holds (e.g. Scudder et al. IQsd ). At higher Mach numbers, where the 
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shock is turbulent, the importance of such electric fields is less clear. Electron E x B drift 
along the shock front will result in periods of energy loss as well as energy gain by th e cross- 
shock potential, and hence no net heating. It has been argued ( IGedalin et al.l 120071 ) that in 



certain cases the shock front may be sufficiently thin (length scales of order c/ujpe) that the 
electrons are effectively demagnetized. One might expect to see electron heating increase 
with shock velocity (or Ma) once this condition becomes satisfied. Examination of solar 
wind shocks suggests that such thin shocks are rare at best and certainly not ubiquitous. 
We see no evidence for an inc rease in electron he ating in SNRs up to shock velocities of 
6000 k m s"^ (IE 0102 2-721 9; Iflughes et aPboOoh . and possibly up to 20,000 km s"^ (SN 
1993 J; iFransson et al.lll996l ). At higher Mach numbers such as those expected in gamma 
ray burst afterglows the convective electron gyroradius may easily reduce to less than the 
electron inertial length, making the cross shock potential a candidate heating electron heating 
mechanism. 



Schwartz et al.l (119881 ) have made a survey of a number of solar wind shock crossings 



observed in situ. They find Te/Tj oc 1/Ma for Ma greater than about 2-3. At lower Ma, there 
is a wide scatter in Tf,/Ti about Tg/Tj ~ 1. At these slower shocks, correlates very well 
with the change in ion velocity squared, suggesting that both are due to the same mechanism, 
presumably the cross shock potential. The switch to Tg/Tj oc 1/Ma at Ma ~ 2 — 3 is possibly 
due to the onset of turbulent shock structure at higher Mach numbers. 

We can explore the conditions required for t he validity of the laminar approximation by 
adopting the criterion of iTidman &: Kralll (Il97ll ) for the existence of a magnetosonic soliton; 



M|/ (M| -l) <Ml< AMI i^s - + 3 + 2 InMs) / (M| - 1 - 2 InM,)^ 



(26) 



The relationship between Ma and Ms predicted by this relation is plotted in Figure 3. The 
criterion above indicates that the laminar approximation breaks down at slig htly lower Mach 



numb ers in the solar wind shocks than indicated by the behavior of Tg/Tj in lSchwartz et al. 
(Il988l ). Magnetic field amplification by about an order of r aagnitude in SNR shocks for 



the 400 km s~^ shocks observed in the Cygnus Loop, where iGhavamian et al.l (120071 ) find 
cor nplete electron- i on eq uilibration, would bring Ma down to the same range as indicated 
by ISchwartz et al.l (119881 ) . possibly suggesting that the cross-shock p otential is at w o rk for 
the lower velocity SNR shocks. At the higher velocity shocks in the ISchwartz et al.l (119881 ) 
sample, which are all perpendicular, an empirical relationship Tg/Tj oc 1/Ma emerges. Such 
a behavior can be consistent with our model if we make the assumption that the cosmic 
rays accompanying solar wind shocks are non-relativistic, suprathermal particles. Then the 
diffusion coefficients take on an extra factor Vg/c, assuming that the cosmic ray velocity is 
proportional to the shock velocity. This extra power of the shock velocity in the diffusion 
coefficient results in Tg oc Vg and Tg/Tj oc l/vg. 
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3.3. Heating versus Damping of Lower Hybrid Waves and the Width of the 

Precursor 



We have discussed whether growth of lower hybrid waves may compete with cosmic 
ray induced magnetic field amplification, but have not yet discussed whether this growth 
rate is sufficient to balance the damping rate of lower hybrid waves by electrons. To an- 
swer this question we compare the electron heating rate from diffusive scattering off the 
lower hybrid waves with the energy input into the lower hybrid waves from the cosmic ray 
turbulence. The electron heating rate per unit area of shock is nefRme>f\\\\tVs/2, where 

fji ^ exp ^— ti;^/2fcyf^g||^ ^ exp (— 2f^mj/f^gme) is the fraction of electrons in resonance 
with the lower hybrid waves (using cj/A; ~ 2vs), is the parallel electron velocity diffu- 
sion coefficient in lower hybrid turbulence, and t is the period of time spent by an elec- 
tron in the turbulence. This time t = l/vg where / is the precursor depth. The electron 
heating is balanced by energy input to the turbulence by cosmic rays with rate 2'yEturbl- 
Putting Eturb = {SE^/Sn) uj^^/fll and >t\\\\ = q^dE'^ky Arn^k^u} we deduce a growth rate 
7 = f Rff B"^ / ^^TT^ifTie^c^ = fn^^/^Q- The kinetic growth rate derived earlier in Equation (7), 
in units of fli, is proportional to ncR/ui. Thus as long as this ratio is comparable to or larger 
than the fraction of electrons that are in resonance with the lower hybrid waves the growth 
rate outlined above will be sufficient to heat the electrons. 

Another constraint on ncR/ui comes from Equation (9). For magnetic field amplifica- 
tion, we require ricRlrii > k\\v\lvslVLi. Since k\\ > l/r^, the gyroradius of cosmic rays at 
injection, ncR/rii > v\/vinjVs ~ 1/1QM\ ~ 10^^ for Ma ~ 10 and Vinj ~ lOf^. Taking 
the maximum growth rate estimated from Equation (9), 7 = MA^incR/2ni ~ Qi/20MA, 
where /j = c/ujpi = va/^i is the ion inertial length, we estimate a characteristic length of 
= 20M^/j ~ 5 X lO^'' cm. This requires a cosmic ray diffus i on co efficient of order 



10 cm s . This is considerably smal ler than th e estimate by iBelll (120041 ). Taking a char- 



acteristic cosmic ray energy of 10^^ eV, iBelll ( 120041 ) finds a typical growth time for magnetic 
field of order 100 years. This would yield a characteristic length scale for magnetic field 
amplification of ~ 10^^ cm for a 3000 km s~^ shock, requiring a cosmic ray diffusion coeffi- 
cient of ~ 3 X 10^^ cm^ s~^. We suspect that our simple estimate reffects the growth rate 
while the shock may be considered quasi-parallel, and that magnetic field amplification slows 
down considerably once it bec omes quasi-perpendicular. Therefore in taking a characteristic 
cosmic ray energy of 10^^ eV, iBelll (120041 ) is taking the lowest energy cosmic rays for which 
the shock may be considered quasi-parallel, and this result may be considered more realistic. 



Further, in iGhavamian et al.l (120071 ) we argued that the depth of cosmic ray precursor 
over which electron heating occurs could not be larger than ~ lO^Vg/rie cm, otherwise neutral 
hydrogen would not survive to encounter the shock front. We suggest here that lower hybrid 
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waves accelerate the small fraction of electrons that happen to be in resonance, and that 
these accelerated electrons communicate their energy to the rest of the thermal population 
by Coulomb collisions, with characteristic time scale 10^° (T/IO^K)^''^ /rig s. Equating this 
to lO^/rie s yields a maximum tem perature of T ~ ^9 ^ .^ (10 ~^)^^^ ~ 5 x 10^ K. This is very 
close to the temperature found in iGhayamian et al.l (120071). 0.3 keV, or 3.5 x 10^ K. Put 
another way, the temperature found in iGhavamian et al.l (120071 ) is consistent with electron 
heating such that neutral hydrogen can survive to encounter the shock front proper. However 
cosmic ray precursors at the small end of the range considered above (~ 10^^ cm) would not 
allow any significant electron coUisional equilibration to occur. Allowing for compressional 
heating of the electrons as they go through the shock, a precursor electron temperature of 
order 10^ K requires a precursor length of ~ lO^fs/ng ~ 10^^ (f^/lOOO km s~^) cm, or a 
minimum cosmic ray diffusion coefficient of D ~ 10^'^ (ti^/lOOO km s~^)^ cm^s~^. 

The electric field in the lower hybrid waves will be given by the limit derived by iKarney 



mm, 



6E = B 



1/3 



UJ 



Ak\c 



B 



1/3 



2k^c 



(27) 



This i s the maximum electric field before ion trapping and heating occurs. iLaming fc Lepri 
(120071 ) and references cited therein demonstrate that when uj/\/2k\\Vte « uj / \/2k±_Vti., ions 
are heated more effectively than electrons above this threshold. In our case, the ions that are 
heated will be the lower energy part of the suprathermal ion distribution refiected from the 
shock, i.e. those below the injection threshold for diffusive shock acceleration in Equation 
(1) or any of its modifications subsequent to the treatment of the reactive lower hybrid wave 
instability in section 2.2. With the wave electric field given by Equation (27), the electron 
momentum diffusion coefficient in lower hybrid turbulence varies as f^, yielding a constant 
degree of heating with shock velocity if the time spent in the turbulence varies as l/ff, which 
would be the case if the cosmic rays are obeying a diffusion law. 



4. Summary 



We have considered in more detail the speculation of IGhavamian et al.l (120071 ) that 
lower hybrid waves generated in a cosmic ray precursor could be responsible for the electron 
heating at collisionless shocks in supernova remnants. We find that there do exist growing 
modes for the resonant or kinetic case, and that the growth rate in this case may be sufficient 
both to survive the damping by electrons and to compete with magnetic field amplification 
by modified Alfven waves. Below a certain Alfven Mach number (roughly estimated to be 
~ 15) the lower hybrid wave growth rate exceeds that of the modified Alfven waves. The 
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modified Alfven wave generation exists for all magnetic field orientations with respect to the 
shock, but is most effective for quasi-parallel case and always generates new perpendicular 
field. Lower hybrid waves, on the other hand, require quasi-perpendicular field geometry 
in order to grow. Thus a schematic picture emerges in which far ahead of the high Mach 
number shock, modified Alfven waves generate perpendicular field, reducing the effective 
Mach number closer to the shock front and thus allowing lower hybrid wave growth to occur 
in a short region before the shock and heat the resonant electrons. A critical Alfven Mach 
number around 15 suggests magnetic field a mplification by about an orde r of magnitude. 



simila r to what a comparison of the surveys of iGhavamian et al.l (120071 ) and lSchwartz et al. 



( 119881 ) would suggest, taking in both cases the shock velocity where (Te/Tp)o ~1 starts to 
break down as that where the laminar shock approximation ceases to hold. 

We have concentrated on the generation of lower hybrid waves, since for these the 
group velocity can be equal to the shock velocity itself, meaning that the waves can stay 
in contact with the shock for long time intervals and in principle grow to large amplitudes. 
However other wave modes that heat electro ns are certainly p ossible, and these, such as the 



Landau damping of kinetic Alfven waves (e .g.lVinas et a.ll2000r). do not require perpendicular 



shocks as lower hybrid waves do. In fact, Bykov fc Uvarov ~ (jl999 ) studied the generic case 



of heating by turbulent modes in the shock precursor and did identify an area of parameter 
space for which a near inverse-square relationship between (Te/Tp)o and shock velocity could 
be accommodated. Our model requires that cosmic ray ions be essentially ubiquitous at SNR 
shocks, with number densities estimated by various means in section 3. In a wider context, 
the idea that cosmic rays are responsible for electron heating at fast shocks reinforces the 
idea that cosmic rays are an intrinsic component of the collisionless shock phenomenon. 
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gram) and by basic research funds of the Office of Naval Research. P.G. acknowledges 
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A. Cosmic Ray Diffusion Coefficients 



The parallel spatial cosmic ray diffusion coefficient is most easily obtained from its 
relation to the pitch angle scattering diffusion coefficient in momentum space . The diffusion 
coefficient in momentum space is expressed most generally as iMelrosd (119861 ) , 



D 



h UJm(X) 



e ■ v(k, p, s)\'^5{ujm - sil 



(27)3 



(Al) 
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where Nn is the number density of wave quanta, Rm is the ratio of electric energy to total 
energy in the wave, such that Rm J Nfj_huMd^k/ (27r)^ = SE'^/Sn, e is the wave polarization 
vector and v is the cosmic ray velocity. For pitch angle scattering by parallel propagating 
Alfven waves, X — a, and so 



AX = h 



—^ + h^ 1 A 



hk\\ 



psina 



With ujm — k\\VA, Rm = {v?\) / {'^c'^) , and e • v = v±/2, 



STT^g^ v\ v"^ sin^ a 



hw 2c2 4 



6{um — sfl — k\iv\\ 



p^sin^a hjjj (27r)3 



p^c^ cos a 



27r 



where we have put, s = 1 and taken ujm << il. 

We now express -Dy in terms of D^a by writing 

/(P, a) = hip) + hip) cos a + ^hip) cos^ a + ... 
and substituting into the diffusion equation 



df ^ df 
+ V 



d 



sin aDn 



dl 
'da 



dt ' ^ dz sin a da 
Upon integrating the result over cos a we obtain, with Vz = vcosa, 

dh , 15/2 , vdh 



dt ^ 3 dt ^ S'dz 



0. 



Multiplying each side by cos a and then integrating over cos a yields 



Idh /^dh v_dh 
3 dt 3 dz ' " 



1^^ - c 

5 dz y_i 



cos a siv? aDf2d{cos a) 



With /o » h » h, h ^ --^ 



3 dz 



J^-^ COS a sin^ aDd{cos a] 



into equation A6 allows the identification 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(AT) 



With equation A3, 



2v^ 
"9~ 



-I -1 



cos a sin^ aD^ad (cos a) 



which when substituted 



(A8) 



(A9) 
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This i s a factor of 2it larger than the equivalent expression given by iBlandford fc Eichler 



( 19871 ). due to a different definition of Um- Where Um oc k.. , Dh oc ^, which evaluates 



to D|| oc vp^^^ or D|| oc vp^^'^ for Kolmogorov or Kraichnan turbulence respectively. If f ~ c, 
the dependence of D\\ on p can usually be neglected. 

The perpendicular spatial cosmic ray diffusion coe fficient has been given in terms of 
D\\ by various authors. Based on numerical experiments, iMarcowith et al.l (120061 ) give = 
?7^"'"'^-D|l where 77 = 6B'^/ (6B '^ + { B)^) and th e cosmic ray distribution function / (p) oc p'^"^. 



Shalchi fc Kourakid (l2007f ) and IZank et al.l (l2006[ ) give oc (^SV^^)^/^ dJ/^ (/azjt; 



\2/3 



from analytic considerations, where I2D is the 2D bendover length scale, the inverse of the 
wavenumber where the inertial range onsets, and consequently has even smaller dependence 
on the cosmic ray momentum than the parallel diffusion coefficient for relativistic cosmic 
rays, and has the same dependence in the nonrelativistic case. 



B. Growth Rate for an Electromagnetic Instability 



For completeness, we give here a treatment of the growth rate due to cosmic rays 
of electromagnetic waves with frequency in the lower hybrid range, and show that it is 
significantly smaller than either the electrostatic instability, or the gr owth of mod ified Alfven 
waves. It is relatively easy to show that the reactive instability of iBelll (120041 ) has higher 
thresholds and lower growth rates as the frequency of the electromagnetic wave increases 
first above the proton gyrofrequency and then above the electron gyrofrequency. Here we 
concentrate on the kinetic instability that might gen erate electr omag netic waves i n the lower 
hybrid range, whistlers, adapting the expressions in iBelll (120041 ) and lAchterberg] (119831 ): 



uov 



+ 47rg 



tOi =F — k 



a 

v±/2 



rip 



JCR ^ ^ 

q k 



(Bl) 



UJ 



(uj-kpw) 



dpi 



Evaluating the two cosmic ray terms 
df 




v_i_— — 27rp±dp±dp^^ 



dp 



[2np±v±f]'^ 



+ k\\v±^^ 2np±dp±dp\\ 



47rp_|_ 



2/2 z 

c k — u . 



7m 



fdp± \ dp\ 



7m 



-ncR (B2) 



and 




df 



{uj — k\\v\\) dp\\ 



2TTp^dp±dp\\ 




-p\/ (7^m^ 



[00 



[1 + ((Pll - mVsY + pl)/i2Kpl)\ 



K+l 




p±{p\\ — mvs)/{'y'^'m?) df 



[UJ 



/cyf ||) 



dpA 



2TTp±dpj_dp\\ 



2np^dp±dp\\ 



(B3) 
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explicitly assuming a k distribution for CR in Equation (B3). After some algebraic manipu- 
lation, an integration by parts, rewriting / as a k distribution and making the substitution 
= P, dP — 2p±dp± (B3) can be written as 



J u-k\\v\\ J \2Kp1 + {p\^-mVsY + P) (mV+pl + P)^ ^ ' 

The integral over dP can be evaluated using a hypergeometic function (Gradshteyn & Ryzhik 
1965; 3.197.1) to give 



ATT{p\\-mVs) {2Kpp'^ B(1,1+k) 

II (m2+p|/c2) (2Kpf+(p,|-mi;.)2)«-i 



/47r(p|| —m 
u)—k\\v\ 



X .F. (2.1;2 + .;l- ^-'^^— -'' jdpii (B6) 

where B(l, 1 + is the beta function. Considering only the kinetic case, using the Landau 
prescription for this integral with the pole at a; — k\\V\\, the imaginary part (i.e. the portion 
relevant for the growth rate) is 

v'j_ (9/ J J \ A- 2 ( I'lTT^ \ i^nPtY ^ 



B(l,l + K 

X 



X 2F1 2,1;2 + ^;1- ^% / ' 27.2 



Setting a; — > a; + i^g in the dispersion relation and taking only the imaginary parts we get 
^ .... + 2.,, ± ^7.^(2.)3,^.„ (W^„ - (,J;S:2).-. B(1, 1 + .)^-l^ 



xF /^2 1-2 + K-1 2^Pt+^X' ^ ^Ci^ (2^-3)r(/.) 
x.F, ^2, 1, 2 + 1 4v/2M3/2p3 r(.-l/2)^^') 



assuming 70;/ /cy -C and including the normalization of / and the factor of A;||/2 that were 
omitted during the evaluation of the integral. Prom this we have 

v. !!£Hgi Mil -""'■) 2K-3 rW p^(2, i; 2 + 1). (B8) 

Hi uj pt + 1) 1 (/t — 1/2) 



Electromagnetic waves in the lower hybrid frequency range are parallel propagating whistlers, 
with 
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and 

dui 2u}/k 

— ~ fBlO) 

dk~ 1 + k^c^/ool^ ^ ' 

We assume duj/dk oc duj/dk\\ ~ Vg and hence uj/k\\ ~ Vs/2 for k « c/cUpg. Thus for «; = 2 
and 

This is significantly smaller than the growth of lower hybrid waves, which is of order 
Qincn/ni. Further, since whistlers carry energy along magnetic field lines, like Alfven waves, 
only for specific shock obliquities will the energy of the waves stay in contact with the shock 
and allow large wave intensities to build up. Electromagnetic waves with frequency above 
the electron gyrofrequency (O and X modes) have phase velocities greater than c, and so 
cannot be excited by kinetic instabilities. 
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Fig. 1. — Schematic illustrations of the amplification of magnetic field by the non-resonant 
modified Alfven waves in the shock precursor in the parallel orientation of the ambient field 
with respect to the shock normal. The evolution of a single field line in an exponential purely 
growing mode is shown. As the field is amplified, the shock becomes quasi-perpendicular and 
the effective Ma decreases, eventually to the point where lower hybrid wave growth takes 
over, allowing a short region of electron heating. 
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Fig. 2. — Same as figure 1 but for an initially perpendicular shock. The evolution of a 
purely growing mode is illustrated. The magnetic field amplification is less strong than in 
the quasi-parallel case, and the shock geometry remains quasi-perpendicular. 
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Fig. 3. — Allowed range of as a function of Ms for the existence of a magnetosonic 
soliton, from Tidman & Krall (1971). The upper hmit is given by the sohd line, the lower 
limit by the dashed line. 



